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^^ I Abstract: The probability of observing a large deviation (LD) in the number of particles 

^^ ■ in a region K in a dilute quantum gas contained in a much larger region V is shown to 

decay as exp[— |A|Ai^], where |A| is the volume of A and AF is the change in the appro- 

r^ ■ priate free energy density, the same as in classical systems. However, in contrast with 

—i I the classical case, where this formula holds at all temperatures and chemical potentials 

*>f-\ ■ our proof is restricted to rarefied gases, both for the typical and observed density, at least 

for Bose or Fermi systems. The case of Boltzmann statistics with a bounded repulsive 

potential can be treated at all temperatures and densities. Fermions on a lattice in any 

dimension, or in the continuum in one dimension, can be treated at all densities and 

^ . temperatures if the interaction is small enough (depending on density and temperature) , 

^ ' provided one assumes periodic boundary conditions. 



^ 



C/3 



^ 



1. Introduction. 



We study the probability distribution of the number of particles in a box A for a quantum 
system in a region V D A described by a grand-canonical ensemble with reciprocal 
C^ , temperature P and chemical potential /i. We are primarily interested in the case where 

the volume of V is much larger than that of A, \V\ » |A|, i.e. in properties remaining 
valid in the thermodynamic limit y — > M with /i, /3 fixed and |A|/|y| — > 0. The particles 
interact with a pair potential v(x — y) so the Hamiltonian of the system is given by 



o, ^ 

^' Hv^^{--A^,-li)+^v{x,-Xj), x,eV (1.1 

i<3 



2 

1=1 

where the Laplace operator is considered with Dirichlet boundary conditions on the walls 



> 

Qs ' of V , which we assume for simplicity to be a cube with side length L (centered at the 

^S] , origin), or with periodic boundary conditions in the Fermionic case at small v considered 

C*"- ' in sec. 7. The analysis in sections 3-6 carries over to the cases of Neumann, periodic 

^-^ , or mixed boundary conditions and for most (reasonable) shapes of the containers V and 

A. This does not imply that effects of the boundary conditions in quantum systems are 
well understood: it just means that in sections 3-6 we do not consider values of /i and 
C^ ■ T where phase transitions are possible. In fact the cases considered in sec. 7. make 

essential use of the periodicity in the boundary conditions. 

The potential v{x) is assumed to vanish for |x| > D and to be bounded, smooth and 
y , stable. The particles will be assumed to obey either Boltzmann, Bose or Fermi statistics. 

Hard core interactions can also be treated, leading to results analogous to those of Sec- 
^ , tions 3-6 but the Fermionic cases in Section 7 would require new ideas in the presence 

of hard cores. For simplicity we do not consider this type of interaction here. 
The number of particles in a region A which, again for the sake of simplicity, we take to 
/K^ ' be a cube of side length £ centered at the origin is given by Na{x) — J^i XA^^i) where 

X — [xi, . . .) denotes a configuration of particles in V and xa is the indicator function 
of the box A {xa{x) = 1 if x S A and xa{x) — otherwise). The dimension of the space 
will be taken d — 3 but our analysis holds in any dimension. The probability of finding 
exactly n particles in the box A is then given by 

n(n) = ^;;';!,,^. (1.2) 

where 5n{x) = 1 if Na{x) = n and otherwise. 
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We note that formally the only difference in 11 (n) between quantum and classical sys- 
tems, is the nature of the probability distribution of configurations x in V . For a classical 
system this is given by a Gibbs measure while for a quantum system it has to be computed 
via the density matrix. 

The pressure P, Helmholtz free energy / and expected density p at temperature /3~^ 
and chemical potential /i are given as usual by 

/3P(/3, y) = lim L-^ log Tr e"''^- , /3/(/3, p) - inf {/3mp - (iP{(i. m)} (1-3) 

L^oo /J 

We shall restrict our analysis to values of /3, /i where the above functions are differentiable 
and the extremum is achieved at a single point p = p{(3, p) so that the expected particle 
density in V 

p=-^P{P,p), ^^±f(p,p) (1.4) 

We call the region in which this differentiability holds the "no phase transitions region" . 
This is a larger region than the region A{G) where the functions P, p, f have been proved 
by Ginibre, see the review [Gi71], to be analytic in z = e^^: a disk in the z-plane with 
radius R{P) > (an estimate for R{P) is quoted later, see (3.4)). The region A{G) is also 
a region where upper and lower bounds on the derivatives ^, -^ can be established: 
Let 

AF{(3, p, po) ''^' /(/?, p) - /(/?, po) - ^4p^ {P Po) 

op p=po , (1.5) 

= [./(/3, P) - P{P, Po)p] - [/(/3, po) - m(/3, Pa)po] 

AF is the difference between the Helmholtz free energy at density p and its linear ex- 
trapolation from pq (which is positive where (1.4) holds if p 7^ po)- We shall prove that 
the probability of finding a density p in A, when the expected density in T^ is po, is given 
asymptotically by 

^g Uin^pe^)^e-''^^^^-P-P''^'" (1.6) 

for large i, provided p— po is small enough. More mathematically this is stated as follows: 

Let /3, p be fixed in the analyticity region A{G) and let po be the corresponding density 
(so that Po = p(/3, p)). If the side I of the box A tends to infinity and, correspondingly, 
the container side also tends to infinity so that L/£ -^ 00 then 

17 limr-^log V n(p^'^)= max -/3AP(/3,p,po) (1.7) 

This holds for all statistics if the interval [a,b] is contained in an interval [— (5o(/3, p), 
(5o(/3,p)] centered at po with 5q{j3,p) > small enough (estimated in (3.4))- 

In Section 5 we obtain (1.7) by proving its finite £ version (which contains various 
finite size corrections). We then prove, in 6, a similar result for Boltzmann Statistics 
with arbitrary (/?, p), (3 > assuming a bounded repulsive (i.e. positive) interaction 
v. this extends the validity of (1.6) far beyond the analyticity region A{G)- We can 
also treat Bose statistics in a region somewhat larger than A{G) but still under the 
very restrictive assumption that e^^'^^^^ < 1 where B is the stability constant of the 
interaction potential v, i.e. X]"<7=i ''^i^i ~ ^j) > —Bn for all (xi, . . . ,x„). Our method 
does not seem extendable to more general Bosonic systems in spite of the strong results 
of Park, sec. 5 of [Pa85]. 

A further extension in Section 7, using completely different techniques (see [GMOO] for a 
review), deals with weakly interacting Fermi systems on a lattice: given /3,po the above 
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theorem holds if the interaction potential is small enough (depending on density and 
temperature) and if periodic boundary conditions are assumed on the container V. The 
results can also be extended to continuum systems in d = 1. They hold for values of /x, 
where the method of Sec. 2-6 fails. 

The above restrictions always exclude the region in /?, /i where genuine quantum phase 
transitions may occur (like superfluidity or superconductivity). Substantial further work 
appears necessary to deal with this regime despite the fact that the ideal Fermi and Bose 
gas (no interactions) can be treated completely by other methods, see [LLS99]. 

2. The Ginibre representation. 

The key technical ingredient in our analysis is the Ginibre representation of the quantum 
partition functions, [Gi71]. We describe here only that part of the formalism which we 
need to derive our results. 

Let J7 — (x, Lu) be a Brownian path starting at a; at t = and returning to a; at a later 
time b: this is a continuous function t -^ uj{t) G M defined for t G [0, b]. The time length 
b will be fixed as 6 = /3 in the case of Boltzmann statistics but it will take the values 
b = P, 2/3, 3/3, ... in the case of Bose or Fermi statistics. 

Functions of 17 will be integrated with a measure dft which is defined as dil — dx-Pxx{duj) 
in terms of the conditional Wiener measures P^ (duj) (see [Gi71] p. 343);''" 



j F{n)dn'il J2 ^ ^^y^ j dx j PiP{du:)F{n) (2.1) 



where, see [Gi71] p. 361, the statistics are distinguished by the upper limit * and the sign 
exponent a as 

{* = 1 Boltzmann statistics 

* = oo, cr = Bose statistics (2-2) 

* = oo, cr = +1 Fermi statistics 

and the integration is over all x G y and a>(i) G V for all t. 

For Ofc = {xk^uik) such that ujk has time length jk[3 we imagine that uok consists of jk 
strings, each of time length /3 with the z-th string denoted by uJk,i{t) = ujk{il3 + t), t G 
[0, /3] . We can consider each bit of string of length /3 as a "particle" and the collection f2 — 
(rii, . . . , r2„) as "trajectory configurations" or as "configurations of particles delocalized 
by quantum indeterminacy" The "energy" of a configuration fl is then defined as 

U{n)^— Yl / viuJk,^{t)-UJk',^'{t))dt (2.3) 

which is consistent with the intuitive delocalization interpretation above in which the 
"number of particles" in the configuration f2 is simply ji + . . . + jn. It is convenient to 
introduce two notions of number of particles of f2 inside A as 



^A(il) = VxA(^fe,^(0)), NA{n)=r^y2 XA{^kAt))dt (2.4) 



t 

(k,i) {kA) 



and note that Nf^{Q^) is an integer while NtviQ.) is generally not, except that Nv{Ql) ~ 
Nv{Ql) is always so. 



This is simply defined to be the measure on the paths starting at x and ending at y in a time 6 formally 
given by Pxy{du>) = Px{duj)S(ijj(b) — y) where Px(du) is the usual Wiener distribution 
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The remarkable representation for the grand canonical partition function (due to Gini- 
bre) is, 

Z(/3,/i;T/) '2/'Tre-''^- = V / ,n ^^m^) ^. (2.5) 

where fl — (r2i, . . . , r2„) and z — e^^ defines the activity. (2.5) makes the quantum 
partition function look like the classical partition function of a gas of closed contours of 
length roughly y/P for all values of z in the Boltzmann statistics and for ze^^^ < 1 in 
the Bosc and Fermi case; i? is a stability constant defined by U{xi, . . . ,a;„) > —Bn,yn. 
Remarkably, the quantum reduced density matrices also admit a natural "classical" 
representation in terms of the above contours but we shall not need such a representation 
here, see [Gi71]. 

3. The cluster expansion for rarefied gases. 

One of the important consequences of the above classical representation of a quantum 
system is that it immediately allows us to take advantage of the techniques developed to 
study Mayer and virial expansions of classical gases. The basic remark is that it is "easy" 
to take the logarithm of a classical partition function. By simple algebraic considerations 
the logarithm is obtained as a formal power series in z and the difficult part is to show 
the convergence of such an expansion. 

The formal power series expression of the grand canonical pressure (3P{(3, /i), defined in 
(1.3) is obtained in the form: 

iogzif3,^i,v) = mp,i^;v) = j2 [ z-'^^m^, (3.1) 

The r.h.s. of (3.1) still looks like a partition function with $"^(il) replacing e~^'^'^ — \ see 
eq. (4.15) in [Gi71]. However, unlike the functions e~^^'' — \ the functions ^^{^) have 
the "cluster property" of "decay" at infinity, i.e. if f2 == (H', fl") consists of trajectories 

n' = {n[,...,n[^,), n" = {n'{,...n':^„) with in',n'Y^\n[,...,n'^,,ni...n:„), 

then 

|$'^(il', fl")| = iid{n',n")>D (3.2) 

where d{A, B) is the distance of two sets A, B and the trajectories are considered here as 
the union of the sets of points they occupy as time varies. This is an elementary property 
that follows from the explicit expression for the functions $"^ in terms of "Mayer graphs" , 
see eq. (4.3) in [Gi71] (see also [GMM71], eq. (4.2) and (4.3) at p. 176, for a similar 
simpler case). 

The functions $"^(^1) are like the functions e"'''^'^ — •• translation invariant in the sense 
that they have the same value for f2 and for a translate of O as long as the two contour 
configurations are inside V: this is a property that will be used to guarantee the existence 
of the limit as L ^ oo. 

Furthermore trajectories that are long give a small contribution to (3.1) at small fugac- 
ities because, see eq. (4.39) in[Gi71], 

/ \dQ\\^^{n,Q_)\<{2-^e-^'"Y-^R-^P+'i-^^ (3.3) 

J N{e_)=p 

Here, q — N{Jl) and w is the maximum of u, \d Q\ means that we use the Bosonic mea- 
sure (which maximizes the integral) whether the system verifies Bose, Fermi or Boltzmann 
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statistics, N{Q) is the number of elementary trajectories composing the configuration 
O and 

R^R{(])^2-^eexp-(3w -2(3B-f3l{(3y'^e''^^C j \v{x)\dx (3.4) 

In (3.4) /(/3) = V2^ is the "thermal length" and C = Ej°li 2"^'.^''/^ see [Gi7I], eq. 
(3.16) (here we made a special choice for the parameter ^ in [Gi71], namely ^ = 2~^e^^"'). 
Since w < oo hard cores cannot be considered. The bound in (3.3) holds uniformly in V 
and it can be considerably improved in the case of Boltzmann statistics, sec [Gi71] eq. 
(3.15). 

Note the asymptotic values of R{(3): R{co) = always and -R(O) = if cf > 2 while 
R{Q) = ro > if d < 2. 

We shall further need the following remark: suppose that the partition function defined 
in (2.5) is altered by inserting a factor f{Q^) which has the property of factoring: 

^(fl', fl") = ¥^(fl')'^(fl") (3-5) 

for all f2', f2" then (remarkably) 



\ogZ,{P,^,,V)'^logY, j ,«^(fl)^e-^^(il)^ 

Consequently if \z\ maxj^ 1(^5(^^)1 < R the series in (3.6) converges uniformly in the size 
L of the container and the sums in the r.h.s. of (3.1) and (3.6) are therefore convergent 
representations of the logarithms of the partition functions (2.5) and (3.6) that define 
them. The bound (3.3) goes back to the theory of the Kirkwood-Salsburg equations and 
of the Mayer and virial expansions, see [Ru69], [Gi71]. 

4. Laplace transform of the probability. 

Turning back to the large deviations problem it is natural to look for properties of 
the Laplace transform (generating function) r(A) of the probability distribution 11 (n), 
defined in (1.2) 

m'^'T.^''-nin)^^-;,^ (4.1) 

n=0 

This admits again a simple Brownian path representation: 

r(A) = ^"=°^^-<^'=" ,^ "' (4.2) 

Note that N\{£1) defined in (2.4) appears here. The similar expression with N\{£1) 
replaced by N\{Q) is a representation of the ratio: 

rpj.g/3A7VA-/3ffv 

4.3 Xr e^f^v ' •' 

which is related to (4.1) as we shall see, but which is interesting in its own right. 
The theory of sections 2 and 3 applies to the expressions (4.1), (4.3): in fact, since 
the functions ip{[l) — e'^'^^^^ — ^ have the factorization property (3.5), eq. (4.1) can be 
written as 
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r(A)=cxp(f: / ^^^^ z-ie^^-^i^^^D^-m^) (4.4) 

ri=0 JVv,(Q) = re 

Let now z be less than i? then {c.f.r. (3.4)), if 

|z||e''^|<i? or X<S{(3,fi)^p-\\ogR-l3^i) (4.5) 

the right side of (4.4) converges to a hmit as V^ ^ oo. This Uniit of r(A) is uniform in V. 
Clearly this is a key point that holds because of the bounds (3.3) and (3.4). Furthermore 
the same convergence bound implies that the argument of the exponential in (4.4) can 
be differentiated term by term with respect to A still yielding convergent series. 

Conclusion: the generating function (4.1) is analytic in A in a small interval |A| < 6((3, ji) 
around 0. In fact one has analyticity in a bigger region of A which contains the infinite 
real semi-axis A < 6{(3,ii). 

We define for each A, |A| = l'^, 

/3P,(/3,M;A)=r'^f: /_^^^ ,«(eW^(il)-l)ci>^(fl)^ 

n=0 «^(^)=„ • 

n=0 NY(0) = n 

The convergence bounds in (3.3) and the translation invariance of the ^^ functions then 
imply that the above expressions have limits as i? ^ oo provided that the size of the 
container V is also such that L/£ -^ oo. (Indeed the integrals over f2 in (4.6) can be 
written as the difference between integrals over completely arbitrary configurations fl 
except for the restriction that they contain one point {e.g. xi) inside A and integrals over 
configurations that intersect the boundary. The latter contribute a quantity of o{£'^) in 
the expression for P and o(l) in the expression for the density).^ The limits of f3Pi and 
pi are given by 



f3{P{f3, p + X)- P{(3, p)) = £ /* ^"(e^''" - 1) <I>^( il) ^ 

p(/3,/i + A) = V I* zV3"n$^(il)^ = 9AP(/3,Ai + A) 



(4.7) 



where the * over the integral means that the point xi = ciJi.i(O) is not integrated over; it 
can be fixed at the origin. 

The approach to the limits is uniform in the parameters (3, ji in any closed region con- 
tained in the analyticity domain A{G). One could also evaluate the finite £ corrections 
and show that if L/£ > 2 (say) then their sizes are of order £~^ and in fact consist of two 
terms of respective orders L^^ and £^^ {i.e. quantities of the order of "surface/volume" 
coming from boundary effects due to the boundaries of V and of A respectively) : however 
one would have to enter into the details of Ginibre's work, so that we just note that (3.3) 
and translation invariance imply that the corrections go to zero as i? ^ oo. 

5. Large deviations in the analyticity region (4.5). 

Given the above information, obtained from [Gi71], the derivation of the large deviations 
results follows the standard path set up in the classical theory of [MS67] . 



In fact a closer analysis would reveal that such terms can be estimated to have size 0(l ^) and 0{l ), 
i.e. that they are boundary terms, see [Gi71] proof of lemma 2.2, p. 366. 
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For (i and ^ — ^q we call, c.f.r. (4.5), pi = Pi{f3, fio', 0) and denote the grand canonical 
averages at chemical potential ^ by (•)^. We estimate the probability 11 that A^a > 
{pi + a)t^, for some small a > 0, by 

n(7VA > [pi + a)l'') < {e^P (NA^(Pe+a)i'))^^ ^ e/3 (p,(0,Po;A)-P,(/3,po;O)-A ip.+a))^^ (5^^) 

for all A > 0. The limit of the coefficient of £'^ in the r.h.s. of (5.1) is, by (4.7), 
(3 {{P{fi, po + X)- P{I3, po) - d^P{fi, po)X - Xa) , which holds for aU A > 0. Minimizing 
over A we get (by definition of free energy), that if a is small enough Il{N/^ > {pi+a)£'^) < 
g-/3a I /2x++o{i ) ^lyQYe x+ is the minimum compressibility {i.e. of the second derivative 
of P with respect to p) in the interval [po — 5, po + S\. This is within the radius of 
convergence of our expansion, provided the value A ^ a/x(0), where the minimum is 

def 

achieved, is such that |A| < S{po,f3), i.e. \a\ < So{P,po) = X+'5(/3, /io) so that the above 
estimates on the Laplace transform hold. 
The case a < can be treated in a similar manner. Hence setting \a\ < 6q{(3,p) we get 

(I) Large deviation property 1: 

U{\Na - pe''] > ai^) < e-/^"^^V2x++o(^^) (52) 

for ae {0,So{(3,pq)). 

The relation (5.2) gives an upper bound on arbitrarily large fluctuations which, since 
the number of possible values of N between —a£'^ and ai"^ is "only" 20^*^ and the total 
probability is 1 implies that inside any density interval of size \a\ < So{(3,po) centered at 
Po there must be at least one value of N whose probability is > 0{£^'^). We write this 
result as 

(II) Large deviation property 2: 

U{Na ^ Poi'^) > 0{a-H-'^) (5.3) 

where N^pl'^ denotes the event \N — pof^] < e£'^. 

This shows that the numbers ^ poi"^ have a "high probability" , and is our second "large 
deviations" result. 

We now consider a density value po + a corresponding to a chemical potential po + Xa, 
with a e {Q,So{po, f3)), and note that the same argument can be applied to the proba- 
bihty distribution generated by replacing e~^^^ in (1.2) by (,-li>-o.Nt,/2^-iiHv ^-liXa.Nf,/2 _ 
Calling n(A'A) this distribution we therefore conclude that 

n(A^A - (po + a^-^) > 0(£-i^-'^)e°(^') (5.4) 

Therefore, with the notation introduced after (5.3), 

n(iVA^(po + a)0 - r^^^-f3H.ef3X^N, 

. ^-KfHp.+ar Tre-P"-e('^^^- ^oie,-) ^ (5.5) 

Tr e-/3ffv 

= m^Hpo + a)£'^)e-^-^(^"+")^^ ^^% ' eOi^''^ 

It e ' ^ 

where we have multiplied and divided by Tre^^^^e^^"^'^ and introduced in the first 

tern in the l.h.s. e^^" ^^ compensating it by Q-^aP(pa+a)e. ^o{(. ) .^^]^j(.]^ jg possible because 

of the restriction A^a ^ (po + fl)^'' in the sum. 
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The first term is between 1 and 0(e ^A '^e '~'^^^ >) by (5.4) and the last ratio can be 
written, again as done above and in the previous section, as 



5.6 "^^ \A^ I ^ '^ ^" ' 'n\ 



dn 



exp(V/ ^^(il) (e/^A.A^ACil) _ 1) lli£) (5.6) 



QJN(n)=n 



n\ 



where the integrals over £2 can be divided into regions with f2 entirely outside the box 
A, which cancel because iVA(il) = 0, and into those partly inside and partly outside, 
which contribute a quantity of o(£''),^ and those entirely inside A, which by the analysis 
in sec. 4 give €^ (i (P(/3, ^^0 + A„) - P(/3, a^o)) = o{(.''). 
Hence we see that 

(III) Large deviations property 3: 

n(iVA^(po + a)0 -e^''^(-^"('"'+")+^('^^'^«+^")-^(^«)) 0{e-H-^e''^''"y) = 

^ g-/3<?''(/(/3.Po+a)-/(po)-9p/(/3,Po)a) Qie^'^ T'^e'^'^^ '■^) 

having used the relations pP{l3, ^iq) — /^MoPo — /3/(/?, Po) if Mo = dp f{l3, po) (see [Ru69], 
eq. (4.17)). 

The smoothness of the functions /, P has been tacitly used in the above discussion. This 
allows us to take advantage of the arbitrariness of e (which we could even choose to be an 
inverse power of i with a small exponent) and to interchange limits and maximizations 
so that (5.7) implies the result stated in sec. 1. Indeed we deduce that in the interval 
[po + a, pq + b]i.'^ there is a value of N whose probability multiplied by g'^-F(/3,po+a)^ jg 
> 0(e(''~")*^'^^ ^(!.~'^) times the exponential in (5.7) and a, b are arbitrary in (0, (5o(/3, Mo))- 
This is the main result claimed in Sec. 1. 

6. Beyond the analyticity region. 

A completely different method can be used to study systems satisfying Boltzmann statis- 
tics at arbitrary chemical potential or Bose gases in the region 6^'^+^''^ < 1. The method 
does not apply to Fermi systems. We illustrate it in the bounded positive interaction 
case: more general cases can presumably be treated along similar lines but we shall not 
attempt to do so here. 

The method is very similar to the one used in classical statistical mechanics by [La72], 
[0188]. This entails comparing r(A) in (4.1) with the partition function Z{j3,p + A; A) 
obtained by replacing Hy in the right side of (4.1) with Hty. More precisely we want to 
prove that \K\^^ logr(A) -^ \K\^^ \ogZ{l3, p + X;A) = P{l3, p + X) which would give us 
directly the desired result. For a classical system this is done by first noting that r(A) is 
(for L — I > D) nothing else than / Z{P, p + X; A\ri)di'{ri). Here 77 is the configuration of 
particles outside A, di'{ri) is the induced Gibbs measure and Z{P, p + X; A\ri) is the par- 
tition function (or normalization factor) for the grand canonical ensemble with chemical 
potential /i + A, in a box A with boundary conditions 77. The proof of the LD formula is 
then basically a version of the proof of equivalence of ensembles, [Ru69]. 

To carry out a similar analysis for the quantum Boltzmann and Bose case it is necessary 
to use trajectories fl instead of configurations. One simply considers probabilities of 
events happening in a fictitious gas of closed paths: the probability of a configuration 
ri = {fli, . . . , ri„) of this gas is, 

z^e-^^<^^Un/Z{p,p,V) (6.1) 



Once more a closer analysis shows that such terms can be estimated to have size 0{£ ^), i.e. that 
they are boundary terms. 
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with the notation of the previous sections. We can now write Q = (f^^, Hext) ^^ 
which the contours in fJ^ are totally or partially inside A and those in [1^^^. do not 
have any points inside A. Then the probability distribution in (6.1) induces conditional 
probabilities Wa(^aI ^ext) verifying the appropriate ^^DLR equations" 

where Z is the appropriate normalization. (Fermions cannot be treated in this way 
because the integration dil is not positive). 

The main estimate on which the analysis will rely is that, in analogy to the classical case, 
the contribution to the partition function from trajectories which cross dK is proportional 
to dK. More precisely there exists a function ??(/3, z) such that 

6 3 ^I^sE/ z'"e-''^(^l^')dil/m!<^(/3,z)<oo (6.3) 

II m J Nv{n)=m. UindK^H) 

for all A C V and for all configurations 17'; the tilde over the integral means that the 
contours in £2 must all intersect A. This is an estimate that can be proved without 
restrictions in the case of a Boltzmann statistics; however in the case of bosons it can 
only be proved under the condition ze'^^^ < 1; for a proof see [Gi71] proof of lemma 2.2, 
p.366.^ 

The probability 11 (iV) that iVA( fl) ~ N can be estimated by considering the probability 
of the same event conditioned on the presence of an external configuration [1^^.^ if the re- 
sulting estimate coincides with (1.6) up to a correction by a factor e*^*- ' independently 
of flej-t the proof of (1.7) will be achieved. 

(IV) Large deviation property 4' 

Suppose that the interaction potential v is positive and bounded. Let /3, /io ^e ii^ the 
region e^^°+^^^ < 1 in the Bose case or arbitrary in the Boltzmann case. We also 
assume that there is no phase transition at /3, /io. Let po be the corresponding density 
(so that po — p{/3,po)), i.e. — i -^ , is continuous at po. If the side £ of the box A 
tends to infinity and, correspondingly, the container side L also tends to infinity so that 
L/C. -^ oo, then 

6 4 limr'^log ^ n(p^'^)= max -/3AF(/3,p,po) (6.4) 

This holds for the Bose system provided the interval [a, h] is contained in an interval 
(0, Pmaa;(/3)) whcrc Pmax{P) is the maximal density that can be achieved as the chem- 
ical potential varies compatibly with the restriction e^^^'^'^^^ < 1. (There is no such 
restriction in the Boltzmann case). 

The assumption of boundedness on the potential does not allow hard cores. However 
the technique below also applies to the hard core case with a tail which is not necessarily 
repulsive: to avoid technicalities we do not discuss this case. Eliminating completely all 



4 



In an unpublished note by Lupini and one of us, (GG), this region was shghtly extended to cover cases 
with ze^ > 1 in systems with hard core interactions. This made use of the fact that hard cores force 
the Brownian paths f2 to be "quite extended" (i.e. 0{a) if a is the core radius) while they have average 

size 0(-v//3): therefore they have a low probability if /3 is small. Hence one can take ze '^ > 1 (by 

0{e°' I ") provided fi is small enough. The analysis could be extended to cover also such cases which, 
however, arc conceptually not really different from the case considered here. 
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6.5 



6.6 



restrictions, i.e. assuming just stability or possibly superstability of the potential, seems 
to require new ideas. 

To prove (6.4) we fix the configuration flg^-j outside A and we look first for a lower 
bound for 11(71). This is obtained by remarking that if £2 C A then for configurations Q 
in which all trajectories are at least a distance D from i9A, \U{il \ He^t)! = 0- 

Hence, if A^ is the region in the interior of A separated from the boundary of 9A by a 
corridor of width D, the numerator of the fraction (obtained from (6.2)) expressing n(n) 
can be bounded below by (/hcAd nA(n)=« ^'''"'"'^"^'"^'^^/"'O which is e^'""" times 
the canonical partition function Q^ for the system with Dirichlet boundary conditions 
(i.e. vanishing boundary conditions) on dAo and no external particles outside the box 
Ad. 

The denominator in (6.2) is bounded above by ZX(^o,/3)e''(^0''')l^^l where Z^ is the 
grand canonical partition function for the system with Dirichlet boundary conditions on 
9 A and no external particles outside the box A; the second factor bounds the contribu- 
tions to the integral from the contours of fl that cross the boundary via the bound (6.3). 
Hence 



Likewise Z^^ gives a lower bound to the denominator of (6.2) while an upper bound on 
the numerator is 

N „ 
fy^ / g/3Mo(W-fc)g-/3;7(0)\g,5(/3,Mo)|aA| (gg^ 

and again we recognize that the integral gives the canonical partition function for N ~ k 
particles in dA with Dirichlet (and no outside particles) boundary condition multiplied 
^^ ^PMN-k)^ Therefore 

Ha. AT < e I I TV max ^ = 

k<N ZlJl3,fio) 

= e''IS'^liV max ^^ ^,,^ ,,,,, e°('^l)(6.7) 

where we have used the fact that the canonical thermodynamic limit (with Dirichlet 
boundary conditions) is uniform in the density in any closed interval contained in (0, Pcp), 
see [Ru69], and furthermore that the maxima in (6.7) are certainly achieved, uniformly 
in the volume, for k > r]\A\ for some 77 > 0. 

In fact the maximum of the r.h.s. of (6.7) is achieved precisely at (5 = by the convexity 
properties of the free energy because 

^ ^ p-f!{fiP,P-S)-f{P,Po)-d,f{l3,po)ip-S-po)) |A| (f. nx 

Where we have used the relation between grand canonical pressure and canonical free 
energy. We have assumed above that there are no phase transitions in the sense that for 
each chemical potential considered there is just one density possible. In particular this 
means that p = — ^ ' and p = '^^-P' and therefore (6.7), (6.8) and (6.5) imply (6.4). 
Existence of phase transitions affect the argument and the results in the same way as 
they do in classical statistical mechanics and we do not discuss the details. 

The above method is purely probabilistic and it relies on two properties: (1) the measure 
dJl is positive and (2) inequality (6.3). Therefore it applies also to the Bose gas case 
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7.1 



7.2 



7.3 



whenever (6.3) can be checked. The inequahty (6.3) does not hold if ze^^^ > 1 so that 
new ideas seem necessary to treat the Bose statistics case in spite of the important results 
of [Pa85] which, by relying on Ginibre's representation of the Gibbs distributions, provide 
a proof of boundary conditions independence. 

7. Beyond the small activity region for weakly interacting fermions on a 
lattice. 

In this section we discuss Fermionic lattice systems. Using the methods developed for 
the study of ground states in Fermionic systems allows one to deal with arbitrary /xq and 
P at small coupling, thus enlarging the region in which LD relations can be obtained, 
at least if one assumes periodic boundary conditions, i.e. regarding V as a torus. (This 
is certainly very restrictive but we note that periodicity of boundary conditions is an 
assumption under which most of the existing results on fermions have been derived.) 

Results can be obtained for all fiQ (i.e. also positive /iq which is not possible with the 
previous method). In the continuum case, however, similar results can be derived only 
in one-dimension or if an ultraviolet cut-off is imposed on k. We shall fix a priori jiq, (3, 
with /3 > 0, and then discuss the large deviations for v small enough (depending on the 
choice of /3, ^o)- 

One can wonder how is it that the low temperature techniques that, even for lattice 
systems, have been proved useful only in 1-dimensional cases (meeting serious, so far 
insurmountable, difficulties in higher dimension) can be of help here in an arbitrary 
dimension. However the low temperature techniques were devised to treat the case 
(3 ~ +00, i.e. ground states. We would encounter the same difficulties in their application 
to large deviations theory if we tried to employ them for the purpose of analyzing LD at 
zero temperature. Applying them at (3^^ > and at weak interaction is, however, not 
difficult at least for lattice systems. 

We consider a d-dimensional square box of even integer side L with periodic boundary 
conditions, and x — (ni, ..., n^), with Ui = — L/2, . . . , L/2. The Hamiltonian is similar 
to (1.1) with the Laplace operator replaced by the discrete Laplace operator. If there is 
no interaction all properties of the system can be obtained through the free Schwinger 
function^ also called propagator, 

where xq £ (-/3,/3) and fco = ^(m+ i), m ^ 0,±1, . . ., and k^ = ^^^ n^ = 0,±1, . . . , 

ztL/2, k — (fci, .., kd), L being the length of the side of V and E{k) — X]i=i(l ^ *^os ki). 
The symbol xq + 0^ means that the value is the limit from the left. 
The "propagator", see for instance Sec. 2 in [BG90] or [BGPS92], is defined in terms 
of the creation and annihilation operators for fermions a^ and of the free Hamiltonian 

T = J2-{^^x^ - Mo) = Jv at (-^^x - Mo)a^ dx as 

g+U,t) =Tr e~(f'-'^^a- e-*^a+ VTr e'^^ 
g_(^,t)=Tre-('3-*)^a+e-*^a,7Tre-^^ ^ ' 

if^ = X — x_' , t>0, which are combined to form a single function: 

Syi^^>-\-g_{-^^-t) if-/3<t<0 ^'■'^> 

And (7.1) is a well known Fourier representation of the function in (7.3), having replaced 
the temperature parameter t with xo for uniformity of notation. Note that the function 
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g{x,t) has to be regarded as a function of t which is defined and periodic, although not 
continuous, in [— /3, /3] (and of x defined and periodic in [—L,L]'^): this is important to 
keep in mind when studying Fourier transforms. 

Note that ki < 27r while fco is unbounded: this unboundedness of the sum is an ultraviolet 
problem "in the temperature direction" and restricting to a lattice has the advantage that 
no ultraviolet problem is present in the "spatial direction" . The ultraviolet problem in 
the temperature direction is somewhat trivial: it is however essential to have no cut-off 
in the ko values. Such cut-offs violate "reflection positivity", i.e. the resulting system is 
no longer Hamiltonian and the physical interpretation becomes unclear. 

We define the "distance on the periodic box V" as dL{xi) = ^sin(^!^) and dp{xo) ~ 
^sin(^^) and dp ^^{x) = {d]:^{xi), ...,di^{xj^},dp{xQ)). This becomes the ordinary dis- 
tance in the limit L ^ oo. 

In the following we shall sometimes write ^ - as J dx. 

For all positive integers a there exists a constant C^, independent of /3, |y |, /ig such that 

where v^^ = 1 + 1/ max{— /io, Z?^^}- 

In fact let h{t) a C°° function which is 1 for i > 2 and for t < 1; we can decompose 
the propagator g into its ultraviolet and infrared parts gmgi] we write 

g{x -y) ^gu{x -y) +gt{x -y) = 

7.5 def 1 ^^hd^^o+0-)+^k^ 1 ^,,_,,,2,,^^!:^^^^f^_ (7-5) 

and we estimate them separately. 
(1) First we check that the bound (7.4) holds for gi{x): if no, ni, ..., n^ are positive 
integers then 

7.6 \{dL{x^r,...dUx^^dp{x„D ■ g,{x)\ < J-^J^Jd^^ ...d^" §,{ k)\ (7.6) 

ko.k 

where dk^, dka denotes the discrete derivative on the lattice of the wave numbers k. Noting 
that the denominator in (7.1) is greater than max{— /LtOj/3~^}j one obtains for gi{x_ — y) 

the bound (7.4) because the sum over fco is finite and that over fc tends to an integral 

over the Brillouin zone which is also finite. 

(2) The next relation that we need is that, for all a there is a constant Ca such that 

C 

\9u{x - y)\ < — — ° -— (7.7) 

l + iy°'\dpj^{x-y)\'^ 

The "only" problem in checking (7.7) is that the sum over fco is unbounded and the 
summand is Oik^ ) so that the sum is improperly convergent if no ~ 0: indeed using 
(7.6) we see that if |no| > 1 the discrete derivative is applied either to the denominator 
(and gives a quantity bounded by (9(fc,^^)) or it is applied to (1 — h{k^)), and dkoh{kQ) 
has compact support. 

Therefore we have to check that the lack of convergence only causes a discontinuity of 
the function g„ at the origin and ±/3 and it does not affect the large distance behavior. 



7.7 
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7.8 



7.10 



7.13 



The non smoothness is just a discontinuity hence it does not affect the size of Qu at the 
origin. A direct check that gu{x) is bounded is obtained by rewriting g„ as 

1 . , pikaxa + k-x( -p(i,\ _ \ i ^ ^ pikoxo 



PV^ {^-^k^ + E{k)-^Ji„){-lko) "■"/? t" -^^0 



fco " (7.8) 



1 d f 



where xC^^o) = 1 if xq > and xl^^o) = otherwise; furthermore ds,y denotes the 
Kronecker delta on the lattice points and we use that (recall that g,gu,9i are periodic, 
although not continuous, functions of Xq in [— /3,/3]). 

..__ ika{xa+0^) i i 

feo = (™+i)27r//3 

The function g is defined by the r.h.s. of (7.8); the first two sums in (7.8) are absolutely 
convergent (the second is a finite sum); hence g„ is finite near xq — 0. Furthermore the 
Fourier transform g( A: ) of p is such that its "Li-norm" {f3V)~^^i^ ^ \'g{k.)\ is uniformly 
bounded in V because (if fco = 'n:l3~^n with n integer) 



\9{k)\< rf^xiM > 1) + ^i ' 1 + T74o,o =» 

^ ^" ^^^ ^ (7.10) 

{tiVr'Y.\9ik)\<c{l + W)P 

ko. k 



for a suitable constant c. Note that instead the Fourier transform of g has a logarithmi- 
cally divergent Li-norm in the above sense because of the last term in (7.8). 

After the above remarks on the nature of the infrared and ultraviolet propagators we 
can go back to the problem of interest to us here. We study the Laplace transforms 
considered in 4, (4.3), and write 

rpj.g-/3(Hv-MoJV)g/3ANA ^ 

= E ^(/3A)" E E XA(27i)...XA(yn)Tr[e-''(^--«^)p(27i)....p(j7„)] ^^'^^^ 

"=0 ' yi Vn 

where p{x) — ipZipZ . Note that since the fermions are on a lattice and in a finite container 
the series over n is a finite sum, by Pauli's principle. 

It is well known that the partition function trace for a Fermionic Hamiltonian can be 
written in terms of Grass-man's variables rj^ = rji ^ , rj^ — Vs x ' ^^ particular, if Hq is 
the kinetic energy (with the discrete Laplacian) 

7.12 Tre-'^C^^-MoAf) J P{dr])e-Uiv) ^ > 

where U{ri) = e J dxdyv{x - y)S{xo - yo)ritv^VyVy , I dx = J^s Iq dxo- The "inte- 
gral" J P(dri)- is defined on monomials of Grassman variables (and extended to general 
functions by linearity) by the anticommutative Wick rule with propagator (7.1). Then 
we can write by (7. 11), (7. 12) 



,PN., JP{dr,)e-^Me^'^^(^^ ^^^^^^ 



' '^° J P{drj)e-uin) 
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7.14 



7.16 



7.17 



7.18 



7.19 



where NAii]) == J2sXa{x)v3mVs,o- 

Remark: had we considered instead of the ratio in (4.3) the quantity 



(7.14) 



where TVa = P~^ J dxg "^g X^ix)'!]'}^ ri~ we would have obtained a different quantity but 
estimates for it can be obtained in a way similar to the ones we explain below: this is 
related to the discussion following (2.4). 

In terms of Grassmanian integrals the logarithm of the Laplace transform of the prob- 
ability (4.1) divided by /3|A| is given by 






7^1og(e^^^^).„-:^log^^^^^^33-,,,;;^ (7.15) 



For a system of fermions on a cubic d-dimensional lattice with Hamiltonian (1.1), in 
which the Laplace operator is replaced by the discrete Laplace operator, then we can 
proceed to developing a cluster expansion at small interaction and proceed to studying 
the large deviations probabilities following the same strategy of Sec. 4, 5 above. 

In fact by the definition of truncated expectation 

log / Pidrj) e^-^^") EE £^{e'^) "il V ^£^(A; n) (7.16) 

where £'^ {A\ n) — £^(Ai, A2, ■ ■ ■ , An)\Ai=A with £{■, •,...) is a suitably defined multilin- 
ear function of its n arguments, we can rewrite (7.15) as 

^ £^{e-^{e^f^'^-l)) (7.17) 



P\A\ 



having exploited that Grassmanian variables anticommutc (so that even monomials in 
Grassmanian fields commute, unlike the Fermionic fields that generate them). More 
explicitly we develop the exponential in powers and obtain 



m\&^^^'^^-^+ f d2LXA{x)vioVs,o-^n) - ^Y.£^{-U;n)] (7.18) 

As usual one can represent the Grassmanian expressions 

U{v) = / dxdyv{x-y)S{xo - yo)vtVxV^Vy_ 

Na{v) = /3A J XA{y)S{y'Hv^ dy 

by the diagrams or ^' graph elements^^ in Fig.l. 



(7.19) 




• 4 
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7.21 



7.22 



Fig.l: The lines "entering (exiting) a point" x represent the Grassrnanian fields 4>'x (respectively ip^ ); 
the wavy line connecting the dots marked x, x' represents v{x — x')S(x'^ — a;"') while the dot in the 
second diagram represents xa(^)<5(^o)- 

We can call the pair P — (x, x') or the single point Q = a^ in Fig.l a cluster. If 

dP = dxdx', dQ = dx, V{P)=v{x-x)S{x° -x'°), 
'•'° Xa(0) - XAix)Six'>), riiP) = r,lr,-^r,l,r,-^,, r,{Q) - r^l^ ^ ^^'^^^ 

and if we set, as above, / dx ■ — J„ dx'^ dx ■ — J„ dx'^ J^s ' ^he (7.18) takes the form 

/m n2 m "2 

{[[V{P,)dP,) ■ (]lxA{Qj)dQ,)-£^{{[[v{m{l[viQ,));ni+n,) 
-i'-^ i=i j=i i=i j=i 

"2>1 

(7.21) 
simply by the multilinearity of the truncated expectations, see [Le87], [BGPS94]. The 
sum over n2 starts from n2 — I because the terms with n2 = cancel because for A = 
(7.17) vanishes. 

We now imagine to join or "contracf ("Wick's contraction") pairs of lines associated 
with clusters (defined above) and with matching orientations in such a way that the 
collection of lines thus obtained, including the wavy lines, form a tree graph connecting 
all points of the clusters. The geometric object so built is called a "spanning tree" and 
contains 2(ni +n2 — 1) solid lines and ni wavy lines. In the graph there will be, therefore, 

de f 

2m = 4ni + 2^2 — 2(ni + ^2 — 1) = 2(ni + 1) lines left "uncontracted" . 

A remarkable algebraic expression for the truncated Fermionic expectations E'^ in (7.21) 
can be developed ( "Lesniewsky's expansion" ) as 

E I'^'Tit) iUgm) ■ (det G^(t)) (7.22) 

T— spanning tree l^T 

where: 

(i) 6 ^ 2; — y if the line I G T joins x and y, 

(ii) t is a family of 2m "interpolation parameters in [0, 1], 

(iii) rT{i)dt is a probability measure on the interpolation parameters whose structure is 

inessential but for the fact that J rT{t)dt = 1, 

(iv) G"^ is a TO X m matrix obtained by considering the m points x^, . . . , x^ ot y , . . . ,y 

into which enter, respectively exit, uncontracted lines and setting G^ ~ tij ■ g{x^ — y .) 

with tij functions of the interpolation parameters such that tiji = u^ ■ u^, for a suitably 

defined unit vectors u^ G i?™. The latter property of the unit vectors is the only property 

that concerns us here. 

The above formula is very convenient because of the Gramm-Hadamard inequality^ 
which applies precisely to matrices that have the form of G^ and bounds their determi- 
nant (for all t) in terms of the rri-th power of the Li-norm of the integrand in (7.1): 

7 23 |detG^(i)|<c5" (7.23) 



This inequality concerns matrices having the form Mj . = («j ■ vj)ff where u^^vj are unit vectors in a 
Hilbert space H and {u ■ v)fj denotes the scalar product in H and it states that the determinant of M is 
bounded by 1 in absolute value. It has, essentially, the simple geometric meaning that the volume of a 
parallelepiped with sides of length 1 cannot exceed 1. An immediate consequence, that we employ here, 
is that matrices Mj = (tij ■ VA)fj (U^ ■ Va)^ with u^, Vj unit vectors in a Hilbert space H and U^, V 
unit vectors in a Hilbert space K also have determinant bounded by 1 (one just applies the previous 
inequality in the Hilbert space H x K). 
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7.24 



7.25 



7.26 



7.27 



7.28 



provided the norm is finite. However, as remarked in connection with (7.9), (7.10) there 
is no cut- off on the momentum feg and the norm diverges ( "logarithmically" ) . Here we 
hit a problem which is serious because apparently the best estimate that we can find 
for this determinant using just the boundedness of the matrix elements is like (7.23) but 
with an extra m!, which would be a disaster. 

Suppose, temporarily, that the propagator g{x — y) is replaced by (/^(x — y) {i.e. we 
impose an ultraviolet cut-off to the model) so that the inequality (7.23) holds: then we 
can quickly develop a proof of a LD result. In fact, assuming (7.23) valid with some 
constant ci, we can bound the ni,n2 term of the series (7.21) by 

-^(/ — 1^ ^U° ( ,|. r+"-'^-(E^("^))'"-^r+^'--^XA(yl)%i) (7.24) 
ni!n2! J v I + i'°'\dpj^{x)\°' ^ 

if we do not perform the integral over the center y of the cluster Qi which certainly 
exists because ^2 > 1: this expresses the cancellation in (7.18) between the terms arising 
from the first truncated expectation when one considers only the contributions from —IJ 
out of the 2" that are generated by expanding the sums in the first expectation (i.e. the 
terms corresponding to n-i =0). 

We see that the difference (7.18) can be given the form 

E X^- fw{X,Y)xA{Y)dXdY (7.25) 

nl,n2>0 '' 

"2>1 

where W{X, Y) is a translation invariant function expressible (as shown implicitly above) 
as a sum of a large number of tree graphs connecting the points in Y U X with X = 
{•S.I, . . . , x^^), Y = {y , . . . ,y ). Furthermore if one fixes the point y € Y , say, the 

integral over the remaining points X, Y^^' = Y/y is bounded, for a suitable choice of C2, 

by 

|iy(X,y^,r(i))|rfXrfy(i) < £"1-14^^+"^ (7.26) 

having used Cayley's tree-counting formula to take into account of the sum over the 
spanning trees and having set e = ^- b(a?)|- We shall show in Appendix Al that (7.26) 
holds as well without imposing an ultraviolet cut-off. 
Hence we can write the difference in (7.18) with the notations of 4, see (4.7), as 

/3(P,(/3,Aio;A)-P,(/3,Mo)=r'^ Y. ^'^ f W{X,Y)xA{Y)dXdY (7.27) 

"2>1 

which has the same structural properties as the first of (4.7) and we can therefore proceed 
in the same way as in 4, 5, with (7.26) playing the role of (3.3), to derive for a system of 
fcrmions on a lattice: 

(V) Large deviation property 5: 

Let (3, /io be fixed arbitrarily and let po be the corresponding density (so that po = 
p{/3,^o)). Suppose that the '^s\v{x)\ — e is small enough, depending on /3,/xo. Lf the 
side £ of the box A tends to infinity and, correspondingly, the container side also tends 
to infinity so that L/£ -^ 00 then 

limr''log y n{pe'^)= max -/3AF(/3,p,po) (7.28) 

<;^oo ^ pela.b] 

pe[a,b\ 
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This holds if the interval [a,b] is contained in an interval [—So(f3, (Iq), (5o(/3, /io)] centered 
at po with (5o(/3,/io) = ^2 v^'^'^'^' > small enough. 

We note that we can replace v{x) by |w(a;)| in the bounds above: hence repulsivity of 
the potential is not necessary and has not been mentioned in the statement above. This 
reflects an essential difi^erence between lattice and continuum systems, see concluding 
remark (7). 

To conclude the proof we refer to Appendix Al where we prove (7.25). The inequality 
(7.26) is studied in [BGPS94]: sec. 3 of that reference is in fact dedicated to a detailed 
analysis of the ultraviolet problem in a case which is more involved than the present one. 
See appendix Al below. 

8. Conclusions. 

(1) In the cases of Boltzmann statistics with arbitrary /?, /io or in the case of Bose 
statistics with ze'^^^ < 1 the Ginibre representation allows us to regard the quantum 
gas as a classical gas of contours, at least for the purposes of computing the partition 
function or the probability distribution of the number of particles in a given region. This 
is so because the integration measure over the trajectories is with a positive definite 
measure dfi. Hence we can apply the arguments used in the general derivation of the 
large deviations formulae from the theory of the equivalence of ensembles of [La72] as 
done in the classical case in [0188]. The applicability of such methods found an early 
application in the work of Ginibre, see [Gi72] p. 362. Naturally we called the results 
"weaker" because we did not get the "total" control on the free energy and the accurate 
treatment of the finite size effects that (3.3) and (3.6) allow us to derive, but the extra 
generality is nevertheless remarkable. 

(2) The above analysis admits a straightforward extension to electron- phonon systems 
because of the remarks in [GGV70], see also [Gi71], p. 420. 

(3) We have also implicitly obtained that the asymptotic behavior of (4.3) (when £ — > oo, 
L ^ Qo and K/l -^ oo) is the same to leading order as that of the r.h.s. of (4.1), in the 
analyticity regions considered. 

(4) The results for Fermi systems, aside from being restricted to lattice systems, con- 
sidered in sec. 7 cover in the fernionic case a rather different region of the (/3, /io)~plane, 
compared to the ones derived in sec. 4. They are valid at any fixed /?, p (or at any j3, /ig) 
provided the coupling is small enough. Thus sec. 7 provides a stronger connection with 
the results for free fcrmions in [LLS99] which are implied by our results in the interval 
of density (— (5(/3, /io), <5(/3, /^o)) in which the latter hold. 

(5) It may seem that the obstacle to deriving a LD relation in the case of the Bose 
or Fermi statistics and in the region where genuinely quantum phase transitions can 
occur is related to the problem of showing boundary condition independence of the main 
extensive thermodynamic quantities. However this problem has been solved in many 
cases: in some cases, indeed, the technical estimates of its solution imply, see [Ro70] and 
[Gi71] p. 365, our equations (3.4) and (6.3), which have been the basis of our results 
in the first six sections of this paper; but in other cases, and we think of the paper 
[Pa85], boundary conditions independence can be proved for very general interactions 
{e.g. bosons with just superstable, if d < 4, or superstable and repulsive interactions, 
if d > 4): but the result does not imply (6.3) and therefore it seems that LD needs 
more than just boundary condition independence. Perhaps (6.3) has to be replaced by 
estimates that involve all particles coherently, as in [Pa85], rather than single ones as 
it is essentially the case in the proofs of (6.3) (which depend on the finite size of the 
individual closed brownian paths). 

(6) The results in sec. 7 are valid for small interactions: how small depends on /3 and 
/io- The only case in which a similar result can be found for a continuum system is the 
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d ~ \ case. If we consider a l-dimensional continuum system we can proceed in the same 
way after replacing sums over x with integrals and using that the kinetic energy grows 
quadratically at oo so that the integral replacing the sum in (7.4) is still convergent: this 
is peculiar to rf = 1 so that the result does not extend to higher dimension unlike the 
lattice case. Of course if an ultraviolet cut-off is imposed on k the above analysis carries 
over to higher dimension: since ultraviolet cut-offs on the spatial momenta k preserve 
reflection positivity the result has some interest because it expresses properties of systems 
which arc still Hamiltonian in spite of a "strange" form of the kinetic energy for large 
P. 

(7) The key bound (7.26) has been rederived here (see Appendix) but it is implicit in 
the analysis in [BGPS94], Sec. 3, where the more complex problem of the ultraviolet 
stability in one dimensional Fermionic systems has been treated in detail. We have not 
simply referred to [BGPS94] for the sake of completeness (the result is somewhat hidden 
there as the concern was about different topics). It would be interesting to prove that 
the inequality (7.23) holds in spite of the fact that the Li norm of the Fourier transform 
of g is infinite; or, alternatively, to show that it does not hold so that the proof that we 
reproduce in Appendix Al for (7.25), (7.26), independent of the validity of the bound 
(7.23), is in a sense optimal. 

Acknowledgements: We are indebted to H.T. Yau and H. Spohn for several discussions 
and examinations of alternative methods. 

Al. Appendix: Proof of (7.25), (7.26). 

We shall represent the Grassmanian fields rj in (7.13) as 77 +V'^ where rj has propagator 
'g{x — y) as defined in (7.8) while V' has propagator, see (7.8), —x{xo)^xQ- 

The Grassmanian integrals in (7.13) will be written as double integrals: for instance the 
numerator in (7.13) will be written 

,-!7(7T+^)+/3A7Va(77+V) 



P{dJq) / P(d^)e-^W+^)+/:'AA/A(„+^) ^^^_^) 

The difference with the corresponding treatment in [BGPS94] is that we use the simplicity 
of the propagator of -0 to perform explicitly the integration over ■0. The integration over 
■0 in (Al.l) can be performed via the formula (7.16) and we replace (7.17) by 

-i- [£l{£l{e-^^^+^^+P^^-^-'+^->)) - £l(£l{e-^''^+^y))] (A1.2) 

The evaluation oi gT(^e~^(ii+i,)+i3XN^(n+i,)^^ jg similar to that oi £T (^e-u(v)+P>^N^(n)-^ -^^ 

(7.18) through (7.19) above. The ST, will be a function expressed as a sum of Grassmanian 
monomials in the fields rj: each such monomial will arise as a "value" of a suitable graph 
according to a (classical) procedure that we proceed to describe. 

Since U and N\ depend upon rj + ip we develop the monomials in U and N\ obtaining 
a sum of 20 monomials in -0 with coefficients that depend on rj (16 from U and 4 from 
A^a): each of such terms can be represented graphically by a graph like those in Fig. 1 
in 7 in which one or more of the solid lines is replaced by a dashed line (hence we get 
16 different graphs from the first and 4 different ones from the second): the dashed lines 
represent the Grassmanian fields -0 while the solid lines represent the fields rj. 

For instance a graph with 3 solid lines and one dashed line will represent TitrJxVx''^x' ^^ 
in Fig. 1 above a dashed line has replaced the solid line entering x' . 

Keeping the notation introduced in 7 we can say that while in the approach attempted 
in §7 wc had just two types of clusters (denoted P or Q) we now have 20 different types of 
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clusters symbolically represented by ^''graphs elements^^ like those in Fig. 1 but with one 
or more solid lines replaced by dashed ones. The different types of graphs elements will 
be labeled with a label j = 1, 2, . . . , 16, 17, . . . , 20: the last 4 will be the ones generated 
by the second graph element in Fig. 1 {i.e. by the graph element with two lines only). 

The truncated expectation fj is computed according to the usual rules: namely we 

must form all possible graphs 7 with n = J2j=i ^j clusters, of which Uj are of type j, 
in which the dashed lines are pairwise contracted or "paired" to form oriented dashed 
lines connecting nodes of graph elements so that no dashed line is left out "unpaired" and 
the graph obtained is a connected graph 7. The wavy lines present in the first 16 graph 
elements must taken into account in checking whether or not a graph is connected. 
The "value" of a graph will be a monomial in the Grassmanian variables 77 rather than 
just a number as it was the case in 7. 

Consider a graph 7 formed from Uj clusters Pji, i ~ 1, . . . ,nj, oi type j = 1, ... ,16 or 
Qj^i, i = 1, . . . , nj of type j = 17, . . ., with a total of 2ndash of dashed lines. Then 7 will 
represent a monomial in the fields 77 that we denote it^{PQ) and which is the product 
of the fields rj associated with the solid lines (all of which remain "unpaired" because 
pairing only concerns the dashed lines, i.e. the fields ip) of the clusters. 

Since the propagator of a pair of fields V';^; V'J is —^x-x'Xi^o ~ ^'o) a factor (— l)"dosh 
has to be added together with all the Kronecker deltas and all the x{' ^ ') functions 
expressing that the various xq — x'q differences must have a definite sign. 

We call A(7) the value of the product of such quantities: it has value 0, ±1 and we 
denote Pj,Qj the collection of the nodes of the clusters Pji or Qj^i respectively. In 
conclusion the value of the graph 7 can be written as 

/ ( 11(11 ^(^^^^))^) • ( n (n^xA(Q,.))^)-.(^Q) A(7) (A1.3) 

•^ j = l i=l J j = 17 i=l J 

The key remark is that the graph 7 cannot contain any closed loop other than the ones 
formed by contracting lines that exit from the same graph element (giving rise to what 
are usually called "tadpoles") because the x functions in the propagator of 4' force the 
Xq 's of the various nodes of the 20 graph types to be in increasing order in the direction 
of the arrows drawn on the dashed lines that are contracted. 

Note in fact that only the graph elements corresponding to the first 16 types contain 

fields computed at exactly equal times (a necessary condition for the possibility of loops) . 

Therefore the number of graphs 7 which are all trees is not bigger than (X]7=i '^j)' 

and we can use (X],=i '^jO'/Ili '^i- < 20^^ ^ in estimating the number of different 
monomials that we can get in the construction. The number of graphs which are not 
trees because of tadpoles is not much bigger (being bounded by 6"n!). 

In this way we express £:T(-g-(7(7j+v.)+/3AArA(77+v.)) = iog£^(-g-c/(7j+^)+/3A7VA(7)+V)) as a 
sum of monomials in the fields that we can ^^clusters^^ and which can be represented 
graphically in a way similar to the one used in Fig.l, Sec. 7, to represent the monomials 
in —U{ri) + (3XN\{ri): namely as in the following figure Fig. 2 

and analytically as 

f ^ 



AlA n p 

■Vx^ ■■■Vx^---Vx[ ■■■Vx'jYldXi)(Y[dz^) 



'=' (A1.4) 



i^l i^l 
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Fig. 2: Illustration of a contribution to 5J in (A1.2) from a cluster generated by collecting all graphs 7 
which have 6 external lines and contain p graph elements of type > 16 with nodes at £1, ■ • ■ , £„• These 
are graphs that contribute with a monomial of the form rj'^ rj~^ jj^ jj~ rj~ jj^ . 



The V function is obtained by collecting and summing contributions from graphs 7 which 
lead to a monomial of the degree 2n in the rj fields and contain p graph elements with 
only two lines. Such contributions will admit, if p > I, the bound 

^1,5 / \^ [SLii ■ ■ ■ 1 3Ln^ 2Li: ■ ■ ■ J 2Lm 2.1: ■ ■ ■ 1 ^p)\d-x_dx_ az_2---dz_p<C2 (Al.b) 

Note that the kernel V may contain some delta functions (or Kronecker deltas) as it 
happens already in the corresponding representation for ~U and X(3N\ described after 
(7.12). 

We do not discuss further details because the technique has been used many times in 
the literature and in particular in [BGPS94]. The validity of the bound in (A1.5) is an 
immediate consequence of the fact that all the graphs that contribute to V are connected. 
This means that in the graph value one can associate a factor v{x — x')Sxo-x' to each 
wavy line that connects x and x' or a factor Ss^s> (while no decay is present in the 
xo — Xg) if the points are linked by a dashed line. The lack of decay in the xq variables is 
not a problem because the xq 's vary in a finite interval [0, /?]. The combinatorics requires 
some care and the analysis is based on the remark following (A1.4) and is also discussed 
in [BGPS94], for instance. 

We can now proceed to the integrations over the rj fields. 

The procedure is the same as the one above. A minor difference is that we have now 
infinitely many graph elements: namely all the ones of the form illustrated in Fig.l and 
in (A1.4). This is not really a big difference because we have the bound (Af .5) which 
will allow us to control the proliferation of graph elements. 

This time we can get, however, graphs with loops and we cannot control their number 
unless we are willing to obtain bounds that grow factorially with the size of the graph. 
To avoid the appearance of factorials we must collect many graphs together; or we must 
compute the expectations by using Lcsniewsky's formula to avoid getting combinatorially 
large contribution. 

The idea is that this time the propagator g of the fields rj, that we have to integrate, 
has a finite Li norm, as pointed out in the comment to (7.10), so that we can apply the 
formula (7.23). 
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Furthermore the second truncated expectation appearing in (A1.2) will not contain any 
term with p > (of course) , but it will otherwise coincide with the expression of the first 
expectation. So that the difference of the two will be given simply by the sum of the 
terms with p > 0. And combining (A1.4) with (7.23) one gets an expression like (7.25) 
and the bound (7.26) with the symbols Xi,... and 2.1, .. . taking the place of the X and 
Y respectively. 

Logically we avoided proving the estimate (7.23) for the full rj integration, and we still 
do not know whether such an estimate is really valid. This has been done by performing 
exactly the (relatively) simple integral over the component ip of the Grassmanian field 
T] = rj + ijj which is responsible for the summability problem described in the comments 
to (7.23). In this way it becomes clear that the logarithmic divergence of the Li norm 
of the Fourier transform of the propagator cannot really cause problems in the estimates 
that we need for the purposes of proving (7.26) and the LD relation 5 of 7. 

It is worth stressing that the problem that we studied in this appendix (originally solved 
in [BGPS94]) is usually considered "irrelevant" in most of the literature: in condensed 
matter Physics it is very common to find that the infrared problems are directly attacked 
by assuming that the propagator contains an ultraviolet cut-off both on the k and the feg 
components.® This is done even in works in which the authors are aware that a cut-off 
on fco is very difficult to interpret because it breaks the reflection positivity of the theory, 
i.e. makes it "non Hamiltonian" . 
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